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We investigate the survival of circularly polarized light in random scattering media. The surprising persistence 
of this form of polarization has a known dependence on the size and refractive index of scattering particles, 
however a general description regarding polydisperse media is lacking. Through analysis of Mie theory, we 
present a means of calculating the magnitude of circular polarization memory in complex media, with total 
generality in the distribution of particle sizes and refractive indices. Quantification of this memory effect enables 
an alternate pathway towards recovering particle size distribution, based on measurements of diffusing circularly 
polarized light. 

PACS numbers: 42.25.Ja, 42.25.Dd. 42.62.Be 


I. INTRODUCTION 

An understanding of scattered radiation is deeply embed¬ 
ded within such scientific disciplines as astronomy, meteorol¬ 
ogy, climatology, and more recently biomedical imaging [1]. 
Measurement of this radiation allows for the remote detec¬ 
tion of object properties, through inspection of both the spec¬ 
tral, and spatial redistribution of energy. Notably, the inclu¬ 
sion of polarization in this analysis serves to enrich the source 
of available information. An area of immense interest con¬ 
cerning such behavior, is the study of optical radiation scat¬ 
tering within biological media [2]. It was first realized that 
through appropriate filtering, polarized light could enable the 
selective imaging of either surface, or subsurface tissue lay¬ 
ers [3]. Since then, similar applications have helped to im¬ 
prove imaging modalities such as optical coherence tomog¬ 
raphy [4], and light scattering spectroscopy [5]. In addition 
to the advancement of existing techniques, the recovery of 
physical properties solely based on the observed transforma¬ 
tion of polarized light is the objective of much research [2]. 
Interestingly, in media consisting of particles comparable to 
the wavelength in size (Mie regime), circular polarization can 
survive many more scattering events than linear polarization 
due to excessive forward scattering. Linear states depolarize 
once the direction of the incident photon stream is eventually 
isotropized [6]. On the other hand, due to the preference of 
shallow angle scattering which generally preserves helicity, 
incident circularly or elliptically polarized states persist on a 
longer length scale. This effect is known as polarization mem¬ 
ory [7-9]. The extent of the polarization memory is strongly 
dependent on the size parameter [10, 11], and on the refrac¬ 
tive index mismatch between the particles and background 


medium [12]. This offers an additional gauge by which to 
characterize medium properties, as measurements of the po¬ 
larization state of multiply scattered light contain this infor¬ 
mation, even beyond the distance required to randomize pho¬ 
ton direction, known as the transport length /,. 


In this article, we address the dynamic relationship be¬ 
tween circular polarization memory and the particles present 
in the scattering medium. Recent studies have shown that 
the survival of linear polarization exceeds that of circular po¬ 
larization in many biological tissues [13-15]. Such behav¬ 
ior is characteristic of smaller Rayleigh scattering particles, 
although these tissues exhibit anisotropic scattering typical 
of the Mie regime. This has suggested an inability to pre¬ 
dict depolarization rates in complex media based purely on 
bulk scattering properties such as the reduced scattering co¬ 
efficient yu', the absorption coefficient p„, and the anisotropy 
g factor [16, 17]. Here, we demonstrate a quantitative de¬ 
scription of circular polarization memory in media consist¬ 
ing of a variety of independent scattering spheres, with total 
generality in particle size, and refractive index distribution. 
Our method uses direct analysis of the single scattering be¬ 
havior attained from Mie theory, and a flexible integral ap¬ 
proach towards characterizing the depolarization rate within 
complex media. This will allow for additional information 
to be gathered from diffusing light measurements, and clarify 
the connection between particle distribution, and the observed 
crossover between scattering regimes. 
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II. HELICITY SURVIVAL PARAMETER 


Mie theory offers an exact solution for scattering of a plane 
wave by a single spherical particle of any size. The relation¬ 
ship between the incoming, and scattered Stokes vectors de¬ 
scribing polarization can be expressed in terms of the scatter¬ 
ing matrix [18]: 
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where the elements of the scattering matrix are dependent on 
the scattering angle 0, and the Stokes vectors are both de¬ 
fined with an orientation aligned to the scattering plane. For 
an incident right-circularly polarized state, with /, = V, = 1, 
Qi = Uj = 0, the Stokes vector for the partial wave scattered 
in the 0 direction, with (ft azimuth is: 
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which is identical for all <p due to the rotational invariance 
of the input Stokes vector. The expression for the degree of 
circular polarization of this partial wave is: 
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while the ensemble average of this quantity, denoted Pc, over 
all scattering angles, weighted by the intensity of each partial 
wave is thus. 
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After any arbitrary number of scattering events n from identi¬ 
cal independent spheres, the value of Pc for all partial waves 
can be shown to follow an exponential decay [8]. Here we 
will define the characteristic helicity survival parameter Nc, 
which is a measure of the number of scattering events required 
to depolarize an incident circular state. The average degree of 
circular polarization as a function of scattering order n is then: 


Pc{n) - e 


-nINc 


( 5 ) 


where Nc can be found from the single scattering behavior 
calculated in Eq. (4) and the natural logarithm of Eq. (5) with 
« = 1. This procedure can be carried out for a monodisperse 
system for any size parameter X - 2nani,lA, and any refrac¬ 
tive index ratio m - njni,, where and n* are the refractive 
indices of the sphere, and background respectively, and a is 
the particle radius. Figure 1 shows the calculated survival pa¬ 
rameter for a range of sphere sizes, and refractive indices. The 
value of Nc can be seen to fluctuate on both large and small 
scales with respect to a change in the size parameter X. 



FIG. 1. (Color online). Helicity survival parameter of perfectly 
monodisperse systems of scattering spheres with size parameter X, 
and refractive index ratio m = 1.15 (dash-dotted black), m = 1.20 
(solid blue), and ra = 1.25 (dashed red). 


III. POLYDISPERSE MEDIA 


For a medium consisting of a range of particle sizes, Eq. (5) 
must be generalized to a distribution of exponentials: 
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where If{X) is appropriately normalized such that, 

IjiX) dX=l . (7) 

I fix) is the fraction of total intensity scattered, on average, by 
particles of size X with corresponding helicity survival param¬ 
eters given by NdX). This function is directly proportional to 
the scattering coefficient of each particle size, I/iX) oc fisiX). 
In other words, the resulting circular depolarization rate of 
a polydisperse system is dependent on not only the helicity 
survival parameters of the constituent particles, but also the 
relative scattering cross sections o-dX), and number densities 
P(X). We then have: 
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and Anally, in terms of the volume fraction distribution fiX): 
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where v(A') is the single particle volume for spheres of size 
parameter X. Equation (6) integrates over the helicity survival 
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FIG. 2. (Color online). Circular polarization memory (Nc/Ni) for a 
range of size parameters, refractive index ratio m = 1.04 (red upper), 
and m = 1.2 (blue lower). Monodisperse media (solid lines), and 
poly disperse media with a 10% Coefficient of Variance in particle 
size (dashed lines). 


parameters of all particles present in the medium, with each 
Nc being calculated from the single scattering behavior. Thus, 
with a database of values stored for A^e(V'), and crj(V'), the 
effective helicity survival parameter for any polydisperse sys¬ 
tem of spheres defined by the distribution function If{X) can 
be found immediately, a key advantage of this approach. 

Narrow polydispersion in scattering media effectively 
washes out the short scale fluctuations of particle size on the 
helicity survival Nc- For example, media with particle num¬ 
ber densities p(X) obeying a normal distribution with a small 
(10%) coefficient of variance in particle size are shown in 
Fig. 2. Here the ratio NdNt is used to show the extent of the 
circular polarization memory, where N, is the number of scat¬ 
tering events required to randomize photon direction i.e. the 
number of scattering events in one transport length I,. Nar¬ 
row polydisperse distributions shown in Fig. 2 are common in 
media such as solutions of polystyrene microspheres in wa¬ 
ter {m - 1.2), where the coefficient of variance is typically 
2 - 10%. These solutions are widely used as simulating me¬ 
dia, or “phantoms”, due to their well defined bulk scattering 
properties. 

Scattering from more complex materials such as biological 
tissue can include contributions from a wide array of particle 
sizes. Disagreement between measurements of circularly po¬ 
larized light scattered from tissue samples, and from matched 
phantoms has been reported [14, 19]. This is likely due in part 
to the presence of both large Mie scattering structures in tis¬ 
sue such as cell nuclei (typically 5-10 jum), as well as smaller 
Rayleigh-Gans scatterers (mitochondira, lysosomes etc. < 0.5 
jum) [20]. These phantom studies often use a monodisperse 


solution with particle size and concentration chosen to match 
the reduced scattering coefficient = //^(l - g), and the scat¬ 
tering anisotropy g, to another more complex medium. How¬ 
ever this is an oversimplification, and fails to reproduce the 
depolarization behavior, as we shall demonstrate. To approxi¬ 
mate scattering from cell nuclei with smaller interspersed cel¬ 
lular organelles, we have calculated Nc for media where the 
volume fraction contributions from small scatterers, and large 
scatterers have been varied. Here, f{X) is chosen to take a bi- 
modal form, with one normally distributed contribution with 
size parameter X - 1.3 + 10%, and another normally dis¬ 
tributed contribution at X = 50 + 10%. The fraction of the 
total volume of scattering material taken by the small scatter¬ 
ers is varied from 0 (nuclei only scattering) to 1 (small or¬ 
ganelle only scattering). The refractive index ratio used is 
m = 1.04 which is typical of biological media with Uj ss 1.42 
(organelles), and n* ^ 1.37 (cytoplasm) [20]. The results are 
shown in Fig. 3. It can be seen that the circular polariza¬ 
tion memory rapidly deteriorates with the addition of small 
scatterers to the medium. For relative volume fractions of 
small particles above 10%, the effectiveness of the medium 
to preserve multiply scattered circularly polarized light is sig¬ 
nificantly diminished. To emphasize the importance of ade¬ 
quately accounting for polydispersion in scattering media, we 
have also calculated the circular polarization memory (Nc/Nt) 
for a set of media with only one dominant size. These points 
on the graph (Fig. 3) represent typical “monodisperse” opti¬ 
cal phantoms with a 10% size distribution about the mean di- 



Relative volume fraction of small scatterers 


FIG. 3. (Color online), (solid line) - Circular polarization mem¬ 
ory (Nc/N,) verses relative volume fraction of small scatterers of size 
X = 1.3 ± 10%, with remaining scattering material composed of 
particles with X = 50 ± 10%. Refractive index ratio is m = 1.04. 
“Monodisperse” phantoms with equivalent macroscopic scattering 
anisotropy, refractive index ratio m = 1.040 (circles), and m = 1.2 
(squares). (Inset) - Matched scattering anisotropy of all three types 
of media. 
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ameter. In each case, these phantoms have been chosen to 
match the scattering anisotropy of the bimodal distribution 
used to model tissue. The circles represent phantoms con¬ 
sisting of particles with m = 1.04. The squares are calcu¬ 
lated with m = 1.2, to represent polystyrene phantoms. The 
branching of these points is due to the multiple possibilities 
in particle size for certain anisotropy g values, and the partial 
absence of data points expresses that the upper limit of scat¬ 
tering anisotropy for polystyrene microspheres is g = 0.93 (in 
water). It is clear from the figure that a significant overes¬ 
timation is made of circular polarization memory when such 
phantoms are used, even if they are chosen such that fig, g, and 
hence /r' are equivalent. The shortcomings of monodisperse 
phantoms evident here are in direct agreement with the previ¬ 
ously mentioned studies [14, 16, 19]. It is important to note 
that the ability of small scatterers to reduce the polarization 
memory of larger constituent particles is heavily dependent 
on the relative scattering cross section cr{X). In the example 
situation shown in Fig. 3, X - 1.3 spheres were chosen be¬ 
cause they are more effective at reducing circular polarization 
memory than smaller (X < 1.1) spheres, which have intrinsi¬ 
cally lower values of Nc- This is due to the trade off between 
faster depolarization, but lower scattering efficiency as size is 
reduced. This demonstrates that depolarization of incident cir¬ 
cular states within polydisperse media is an intricate process 
dependent on the specific size and refractive indices of each 
contributing particle. Equation (6) can be generalized further 
to include variations in particle refractive index mismatch m: 
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where IfiX, m) now represents the fraction of total intensity 
scattered by particles for which Nc - NdX, m). The effective 
helicity survival parameter of the polydisperse medium can 
once again be found by inverting Eq. (5) after evaluating this 
integral. To further highlight the flexibility of the approach 
presented here, Eig. 4 shows the magnitude of circular polar¬ 
ization memory for a wide range of size parameters X, and 
refractive index ratios m. Each point on the image represents 
NcINt for a 10% particle size distribution about the mean size 
parameter X. Incident circularly polarized light is seen to be 
preserved most in media where both \m - 11 << 1, and particle 
size is large. Such properties are consistent with large cellular 
structures in tissue. A complex cellular model could be pro¬ 
posed by the integration over all types of organelles present, 
including vacuoles with m < 1. However, the closely packed 
nature of scattering centers in many biological tissues may 
cause discrepancies from such models when dependent scat¬ 
tering arises, as shown in Ref. [21]. Thus, care is required 
when generalizing these findings which are a result of the in¬ 
dependent scattering approximation for spherical particles. 
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FIG. 4. (Color online). Contour plot of circular polarization memory 
(Nc/N,) for size parameter 0.1 < A < 80, and refractive index ratio 
0.85 < m < 1.5. Each point represents the circular polarization 
memory for a narrow poly dispersion (10% coefficient of variance) 
about the mean size parameter X. 


IV. SUMMARY 


In summary, we have presented here a method for charac¬ 
terizing circular polarization memory in polydisperse media, 
where the particle distribution has been demonstrated to intri¬ 
cately determine the decay rate of incident circularly polarized 
light: Narrow polydispersion about one dominant size param¬ 
eter has the effect of smoothing out the short scale fluctuations 
seen with theoretical solutions containing perfectly identical 
centers. This type of polydispersion is typical of phantom me¬ 
dia such as polystyrene microspheres, and demonstrates that 
models using pure monodispersion can differ somewhat from 
realistic media where statistical variations in particle size are 
unavoidable (see Eig. 2). In complex media, where size distri¬ 
bution spans from the Rayleigh regime to the Mie regime, the 
extent of circular polarization memory is shown to be greatly 
impacted by only a small volume fraction of sub-wavelength 
scatterers. In our particular example (Eig. 3) we showed that 
circular polarization decays within approximately two trans¬ 
port lengths {Nc/Nt < 2) for media consisting of over 20% rel¬ 
ative volume fraction of small 1.1 < A < 1.5 particles. With 
the method put forward here, it is now possible to rapidly cal¬ 
culate the helicity survival parameter Nc for a medium with 
any general form of polydispersion in both size parameter X, 
and refractive index ratio m. Thus, recovery of this parame¬ 
ter from degree of polarization measurements of semi-diffuse 
light can provide an alternate means towards characterizing 
particle distribution. Eurthermore, although circular polariza¬ 
tion memory is a manifestation of anisotropic scattering, we 
have demonstrated the important distinction that it cannot be 
characterized solely by the scattering anisotropy g parameter. 









5 


We acknowledge support from the William Evans fellow¬ 
ship provided by the University of Otago for Professor Steven 
Jacques, The Dodd-Walls Center for Photonic and Quantum 
Technologies, and the Department of Physics at the Univer¬ 
sity of Otago. 


* http://www.biophotonics.ac.nz 

[1] D. A. Boas, C. Pitris, and N. Ramanujam, Handbook of 
Biomedical Optics (CRC, Boca Raton, FL, 2012). 

[2] N. Ghosh, and I. A. Vitkin, “Tissue polarimetry: concepts, chal¬ 
lenges, applications and outlook”, J. Biomed. Opt. 16 , 110801 
( 2011 ). 

[3] J. M. Schmitt, A. H. Gandjbakhche, and R. F. Bonner, “Use of 
polarized light to discriminate short-path photons in a multiply 
scattering medium”, Appl. Opt. 31 , 6535 (1992). 

[4] J. F. De Boer, et ai, “Two-dimensional birefringence imaging 
in biological tissue by polarization-sensitive optical coherence 
tomography”. Opt. Lett. 22, 934 (1997). 

[5] R. S. Gurjar, et al, “Imaging human epithelial properties with 
polarized light-scattering spectroscopy”, Nat. Med. (7), 1245 
( 2001 ). 

[6] M. Xu, and R. R. Alfano, “Random walk of polarized light in 
turbid media”, Phys. Rev. Lett. 95, 213901 (2005). 

[7] F. C. MacKintosh, et al, “Polarization memory of multiply 
scattered light”, Phys. Rev. B 40, 9342 (1989). 

[8] M. Xu, and R. R. Alfano, “Circular polarization memory of 
light”, Phys. Rev. E 72, 065601 (2005). 

[9] Y. L. Kim, et al, “Circular polarization memory effect in low- 
coherence enhanced backscattering of light”, Opt. Lett. 31, 
2744 (2006). 


[10] D. Bicout, et al., “Depolarization of multiply scattered waves 
by spherical diffusers: influence of the size parameter”, Phys. 
Rev. E 49 , 1767 (1994). 

[11] W.Cai,efa/., “Analytical cumulant solution of the vector radia¬ 
tive transfer equation investigates backscattering of circularly 
polarized light from turbid media”, Phys. Rev. E 74, 056605 
(2006). 

[12] A. D. Kim, and M. Moscoso, “Influence of the relative refrac¬ 
tive index on the depolarization of multiply scattered waves”, 
Phys. Rev. E. 64 , 026612 (2001). 

[13] S. Alali, et al., “Quantitative correlation between light depo¬ 
larization and transport albedo of various porcine tissues”, J. 
Biomed. Opt. 17 , 045004 (2012). 

[14] M. Ahmad, etal., “Do different turbid media with matched bulk 
optical properties also exhibit similar polarization properties?”, 
Biomed. Opt. Express 2, 3248 (2011). 

[15] V. Sankaran, J. T.Walsh, and D. J. Maitland, “Comparative 
study of polarized light propagation in biologic tissues”, J. 
Biomed. Opt. 7, 300 (2002). 

[16] M. R. Antonelli, et al, “Impact of model parameters on Monte 
Carlo simulations of backscattering Mueller matrix images of 
colon tissue”, Biomed. Opt. Express 2, 1836 (2011). 

[17] N. Ghosh, et al., “Anomalous behavior of depolarization of 
light in a turbid medium”, Phys. Lett. A 354 , 236 (2006). 

[18] C. E. Bohren, D. R. Huffman, Absorption and Scattering of 
Light by Small Particles (Wiley, New York, 2008). 

[19] V. Sankaran, et al., “Comparison of polarized-light propagation 
in biological tissue and phantoms”. Opt. Lett. 24 , 1044 (1999). 

[20] J. M. Schmitt, and G. Kumar, “Optical scattering properties 
of soft tissue: a discrete particle model”, Appl. Opt. 37 , 2788 
(1998). 

[21] V. Sankaran, J. T.Walsh, and D. J. Maitland, “Polarized 
light propagation through tissue phantoms containing densely 
packed scatterers”. Opt. Lett. 25 , 239 (2000). 



